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1.  INTRODUCTION 


The  research  described  in  this  report  was  carried  out  under  Grant  N00014-89-J-3104,  part 
of  the  University  Research  Initiative  (uri5002)  in  the  time  period  15  July  1989  to  14  July  1992. 
It  had  as  its  goal  the  investigation  of  a  general  class  of  thermoacoustic  devices,  of  arbitrary 
phase  in  the  impedance  of  the  sound  field.  This  research  was  to  understand  the  sensitivities  of 
thermoacoustic  devices  to  this  phase,  to  try  to  optimize  it,  and  start  to  experiment  with  practical 
devices  with  optimized  phase.  Along  these  lines,  we  also  experimented  with  resonators  capable 
of  efficiently  applying  such  an  arbitrary  phase  to  a  thermoacoustic  structure.  Because  this  grant 
was  part  of  the  University  Research  Initiative  and  as  such  designed  to  strengthen  the  teaching 
and  research  infrastructure,  Section  3  of  this  report  summarizes  the  infrastructure  development 
resulting  from  this  grant  at  George  Mason  University. 


The  following  persons  participated  in  the  research: 


Graduate  students 

A.  Koren,  M.S.  student  in  Applied  Physics,  graduated  Spring  1992,  will  continue  in  the 
new  Ph.D.  program. 

M.  Hutchison,  M.S.  student  in  Applied  Physics,  has  thesis  yet  to  finish. 

J.  Velazco*,  M.S.  and  Ph.D.  student  in  Electrical  and  Computer  Engineering. 


Senior  personnel 

P.  H.  Ceperley,  Principle  Investigator. 


‘Mr.  Velazco  received  no  direct  support  from  the  grant.  He  worked  with  the  P.I.  on  a  spin  off 
of  the  rotating  waves  research,  a  study  of  rotating  waves  in  electromagnetic  applications. 
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2.  RESEARCH  PROJECTS 


The  following  research  projects  were  undertaken  during  the  report  period: 

A.  General  thermoacoustics: 

1 .  Development  of  a  thermoacoustic  simulator  apparatus 

2.  Computer  study  of  the  effect  of  phase  on  thermoacoustic  devices 

3.  Measurements  of  properties  of  various  thermoacoustic  structures 

4.  Development  of  a  simple  circuit  model  for  thermoacoustic  devices 

B.  Experiments  on  rotating  waves 

C.  Development  of  a  thermoacoustic  cryocooler 


2.1  Development  of  a  Thermoacoustic  Simulator  Apparatus 


Figure  1  snows  the  three  thermoacoustic  simulators  built  and  tested.  These  devices 
generated  a  sound  field  adjustable  both  in  phase  and  amplitude  for  use  in  testing  thermoacoustic, 
porous  structures.  The  first  device1,2,  labeled  A,  was  operated  in  the  0.3  to  70  Hz  range  and 
produced  pressure  amplitudes  up  to  200  Pa.  Its  design  also  included  hot  and  cold  heat 
exchangers  to  impress  a  temperature  gradient  across  the  thermoacoustic  structure  and  the 
pressure  transducers  which  were  connected  to  a  Keithley  570  data  acquisition  board  and  a  386 
PC  type  computer.  After  preliminary  processing,  the  resulting  data  was  displayed  in  MathCad 
for  analysis.  This  simulator  was  used  both  to  measure  thermoacoustic  properties  of  various 
porous  structures  and  to  provide  the  acoustic  and  thermal  environment  to  simulate  the  operation 
of  various  thermoacoustic  devices.  I  will  discuss  the  results  of  such  measurements  in  Section 
2.3  and  such  a  simulation  in  Section  2.6. 
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Simulation  devices  B  and  C  were  built  to  provide  increased  pressure  amplitude, 
up  to  1 1  kPa  rms,  primarily  for  use  in  simulation  and  development  of  a  thermoacoustic 
cryocooler.  Devices  B  and  C  also  have  the  advantage  of  being  an  acoustically  cleaner 
system,  producing  pure  plugs  of  moving  gas,  while  adding  the  potential  problems  of 
leaking  and  friction  associated  with  pistons. 


speaker 


^^heat  exchanger 
porous  structure 


piston 


- —  speaker 


piston 

pressure  transducer 
thermocouple 

layered  mylar  structure 
—  insulation 


B 


Fig.  1.  Thermoacoustic  simulators. 
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2.2  Computer  Study  of  the  Effects  of  Phase  on  Thermoacoustic  Devices 


The  bulk  of  thermoacoustic  devices  previously  studied  use  standing  wave  acoustic  fields. 
Standing  waves  have  a  ninety  degree  phase  angle  between  the  pressure  and  velocity  amplitude 
in  the  thermoacoustically  active  structure.  The  author  of  this  report  has  long  held  that  traveling 
acoustic  waves,  with  approximately  a  zero  degree  phase  angle,  would  be  superior  for 
thermoacoustic  applications.  Under  this  grant,  a  computer  study  has  been  made  comparing  the 
performance  of  several  thermoacoustic  devices  as  a  function  of  phase.  While  this  calculation 
uses  equations  that  have  been  previously  derived3,  these  equations  are  quite  complex  and  beg  for 
a  careful  interpretation.  Figure  2  shows  a  typical  result4  of  our  studies. 
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Fig.  2.  Typical  output  of  a  computer  study  on  thermoacoustic  devices.  The  left  one 
shows  the  lines  of  constant  gain  and  constant  Carnot  efficiency  on  a  grid  of  channel 
width  (distance  between  layers  in  a  layer  structure)  versus  temperature  gradient  for  a 
traveling  wave  heat  engine.  The  right  one  shows  the  Carnot  efficiency  and  power  gain 
versus  phase  for  various  channel  widths.  Because  the  right  graph  has  acoustical  phase 
as  a  variable,  it  includes  the  whole  range  of  traveling  and  standing  wave  heat  engines. 


2.3  Measurements  of  Properties  of  Thermoacoustic  Structures 


There  are  other  good  methods  of  measuring  physical  and  acoustic  properties  of  materials 
such  as  those  discussed  by  the  authors  of  References  5  and  6,  however  the  thermoacoustic 
simulators  are  able  to  measure  these  properties  in  the  same  acoustic  and  thermal  environment 
as  is  encountered  in  an  actual  thermoacoustic  device,  and  are  able  to  provide  a  completely 
adjustable  acoustical  field.  As  we  reported  in  Ref.  7,  we  used  the  first  thermoacoustic  simulator 
to  measure  the  thermoacoustic  properties  of  the  following  thermoacoustic  structures: 

►  Various  grades  of  steel  wool,  parallel  and  perpendicular  orientations, 

►  One  grade  of  filter  sand, 

►  One  structure  made  of  stainless  steel  shim  stock  and  stainless  steel  spacer  wires,  and 

►  Several  structures  made  of  spaced  mylar  sheets. 


The  left  side  of  Fig.  3  shows  an  end  view  of  one  of  the  porous  thermoacoustic  structures  that 
was  fabricated  and  tested.  The  right  side  of  Fig.  3  shows  a  graph  of  the  thermal  roll  off 
frequency  (the  inverse  of  2x  times  the  equilibrating  time  constant)  plotted  against  the  flow 
resistance  of  these  structures.  The  line  indicates  the  theoretical  ratio  for  parallel  plate  structures. 

More  work  is  yet  to  be  done  to  understand  the  discrepancy  between  the  theoretical  line  and  the 
data  points  in  Fig.  3. 


Performance  of  various  porous  media 

one  atm.  air 


Fig.  3  Testing  of  thermoacoustic  structures:  photo  of  the  end  of  a  layered  structure  and 
thermal  roll  off  frequency  plotted  against  flow  resistance  for  various  structures  measured. 


2.4  Development  of  a  Simple  Circuit  Model 


The  design  of  practical  devices  is  usually  expedited  by  the  use  of  simple  models 
which  give  clear  understanding  of  the  design  tradeoffs.  Figure  4  show  a  simple  electrical 
circuit  model  which  we  have  developed  under  this  grant.  This  circuit  models  the 
interaction  of  acoustical  and  thermal  energies  of  a  small  fluid  parcel  in  a  thermoacoustic 
device.  It  also  is  useful  as  an  approximation  for  whole  thermoacoutic  device,  providing 
parameters  such  as  the  absolute  temperature  and  the  thermal  resistance  do  not  vary 
excessively  in  the  device.  The  equivalent  or  analogous  quantities  are  shown  in  Table  1. 


Fig.  4.  Simple  electric  circuit  model  for  thermoacoustic  processes. 
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Available  power  input  (MW/m  ~  3) 


TABLE  1. 

Equivalent  parameters  between  electrical  circuit  model  and 
thermoacoustics 


Electric  Circuit 

Parameters 

Equivalent  Thermoecouetic 

Per ametere  (per  m3) 

electrical  current 

flow  of  available  work.  (W/m* 

voltage 

abeolute  temperature  «) 

capacitance 

conet  ant  pr  a  a  aura  heat  capacity 

per  unit  volume  U/X-m3] 

reeietance 

thermal  reeietance  (W/K-m*) 

alact.  power  generated 

acouetical  power  abeorbed 

by  voltage  eourca 

(or  ganer  at  ad)  (W/m* 

elect,  power  generated 

thermal  available  energy  deneity 

by  current  aource 

expended*  (W/m*) 

‘The  thermal  current  deneity  (averaged  over  an  ecouatic 

cycle)  flowing  from  hot  to  cold  (urate  -W/m*) 

equale  the  thermal  available  energy  deneity  expended  multiplied  by  the  abeolute  temperature  and  divided 
by  the  temperature  gradient  (dT/dx). 

We  reported  in  Ref.  8  a  comparison  of  results  obtained  from  this  simple  model 
with  those  obtained  from  an  exact  computer  simulation  reported  in  Section  2.2.  As 
shown  in  Fig.  5,  agreement  is  very  good  for  thermal  time  constants  (for  gas  in  the 
porous  structure  to  equilibrate  with  the  solid)  less  than  and  roughly  equal  to  an  acoustic 
period.  Since  traveling  wave  devices  are  best  designed  with  short  time  constants,  they 
are  well  modeled  by  this  simple  circuit.  Standing  wave  devices  generally  have  longer 
time  constants  and  consequendy  are  not  as  well  modeled  by  this  circuit.  Depending  on 
the  exact  parameters,  standing  wave  devices  may  require  an  additional  component  or  two 
in  the  model. 
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Fig.  5.  Comparison  of  results  from  the  simple  electric  circuit  model  with  those  from  a 
computer  simulation. 


2.5  Experiments  on  Rotating  Waves 


The  thermoacoustic  simulators  discussed  above  are  convenient  for  providing  an 
acoustic  field  for  thermoacoustic  research,  however  due  to  the  low  efficiency  of  speakers, 
the  simulator  would  be  impractical  for  use  in  the  operating  thermoacoustic  device.  In 
that  case,  the  acoustic  field  would  most  likely  be  provided  by  a  suitable  resonator.  A 
regular  resonator  will  provide  waves  with  a  90  degree  phase  angle  in  the  impedance. 
Rotating  waves  are  required  to  provide  zero  phase  angle  and  if  mixed  with  standing 
waves  can  alternatively  provide  an  arbitrary  phase  angle13.  Rotating  modes  can  be 
created  by  the  proper  combination  of  degenerate  or  nearly  degenerate  standing  wave 
modes.  Our  group  has  been  studying  rotating  waves,  both  as  a  general  wave  phenomena 
and  as  applied  to  thermoacoustic  devices.  While  in  some  sense  rotating  waves  have  been 
understood  for  years,  the  knowledge  is  not  widespread  and  is  usually  only  mathematical. 
We  feel  that  rotating  waves  are  best  presented  in  a  very  mechanistic  or  visual  way  as 
described  in  Refs.  9  and  10.  We  have  also  been  experimenting  with  rotating  waves,  in 
the  two  resonators11  shown  in  Fig.  6.  The  first,  simple  cylindrical  resonator  labelled  A, 
is  of  greater  interest  in  general  physical  acoustics,  while  the  second  labelled  B  is  of  most 
interest  in  thermoacoustic  devices.  In  both,  we  have  observed  rotating  modes. 

We  have  also  experimented  and  achieved  single  point  excitation13'15  as  required 
for  thermoacoustic  traveling  wave  devices.  Because  a  rotating  wave  is  the  superposition 
of  two  independent  standing  wave  modes,  normally  two  exciting  sources  are  required  to 
create  a  rotating  mode.  However  with  the  proper  perturbation  in  the  resonator,  a  single 
source  can  suffice11.  We  have  also  observed  the  rotating  wave  shock  formation  shown 
in  Fig.  7  at  amplitudes  greater  that  1  kPa,  which  is  an  interesting  physical  acoustics 
phenomenon  in  itself.  Spiral  wave  fields  were  also  generated  and  mapped  out12.  These 
are  of  general  interest  because  they  do  not  propagate  along  a  straight  line  from  source 
to  receiver,  even  in  an  isotropic  medium.  They  are  the  wave  type  emitted  by  propellers. 

Also  as  a  spin-off  of  this  research,  is  the  work  by  J.  Velazco  and  the  author  in 
the  area  of  rotating  waves  in  electromagnetic  microwave  cavities  and  their  applications16. 
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acoustical  driver 

cylindrical  resonator 


pressure  transducer 


Concentric  Ring  Resonator 


Fig.  6.  Resonators  used  in  the  rotating  wave  experiments.  The  right  one  is  a  concentric 
ring  resonator,  while  the  left  one  is  a  simple  cylindrical  resonator. 


Fig.  7.  Pressure  versus  time  in  a  rotating  wave  experiment  showing  steeping  or  shocking 
of  the  wave  fronts. 
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2.6  Development  of  a  Thermoacoustic  Cryocooler 


In  order  to  demonstrate  the  utility  of  understanding  phase  in  thermoacoustics,  we 
have  undertaken  to  develop  a  thermoacoustic  cryocooler  based  on  our  new  knowledge. 
This  has  involved  the  test  setups  or  simulators  shown  in  Fig.  1,  labelled  B  and  C,  as  well 
as  some  of  the  thermoacoustic  structures  discussed  in  2.3  above  fabricated  from  mylar. 
The  thermoacoustic  structure  in  the  simulator  is  shown  in  Fig.  8.  While  we  observed 
cooling  in  the  5  to  50  degree  C.  range,  this  was  considerably  less  than  expected. 
Acoustical  streaming  due  to  vibration  of  the  regenerator  material  is  suspected  and  will 
be  further  investigated. 


Fig.  8.  Thermoacoustic  cooler  structure  in  simulator,  as  tested. 
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3.  CONTRIBUTIONS  TO  THE  INFRASTRUCTURE  AT  GEORGE  MASON 
UNIVERSITY 


One  of  the  benefits  of  this  research  for  both  the  DoD  and  the  public  is  that  it  has 
materially  reinforced  the  teaching  and  research  infrastructure  of  a  developing  university. 
George  Mason  University  was  started  only  a  few  decades  ago  and  is  one  of  the  few 
relatively  young  state  universities  in  our  country.  It  has  only  recently  been  developing 
graduate  programs.  This  grant  has  assisted  that  development  by  paying  for  the 
development  of  a  physical  acoustics  laboratory  associated  with  both  the  Department  of 
Electrical  and  Computer  Engineering  and  the  Physics  Department.  It  furthermore  has 
supported  the  principle  investigator  and  his  students  to  do  research  in  physical  acoustics. 
It  has  also,  in  part,  helped  to  make  possible  the  launching  of  a  new  Ph.D.  program  at 
this  school  in  the  area  of  Computational  Science  with  a  strong  emphasis  in  Computational 
Physics.  Because  this  university  is  located  Fairfax  County,  Virginia  (a  major  part  of  the 
Washington  D.C.  suburbs)  and  therefore  is  very  close  to  many  Naval  and  DoD  facilities, 
these  new  educational  and  research  facilities  are  available  to  and  will  serve  Naval  and 
DoD  personnel.  The  principle  investigator,  as  well  as  the  other  faculty  have  students 
who  are  employed  by  the  Navy  (at  NRL  and  David  Taylor,  for  example)  as  well  as  by 
the  other  armed  forces. 

Reference  17  discusses  color  visualization  of  wave  and  fields  surfaces  done  at 
George  Mason  University.  The  grant  made  possible  the  purchase  of  the  Hewlett  Packard 
color  desk  jet  printer  used  for  this  work. 
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4.  SUMMARY 


Research  during  the  time  15  July  1989  to  14  July  1992  has  been  directed  towards 
developing  improved  thermoacoustic  heat  engines,  utilizing  optimal  phasing  in  the 
impedance  of  the  sound  field.  More  specifically,  thermoacoustic  simulators  were 
designed  and  built  to  allow  testing  of  thermoacoustic  structures.  A  computer  simulation 
of  thermoacoustic  devices  was  made  to  study  the  sensitivity  to  phase  in  several 
representative  thermoacoustic  devices.  The  thermoacoustic  properties  of  a  number 
porous  structures  were  measured,  including  ones  made  of  steel  wool,  sand,  layered 
stainless  steel,  and  layered  mylar.  A  simple  equivalent  circuit  for  thermoacoustic  devices 
was  developed  and  verified  with  computer  simulation.  Rotating  waves  were  explored  as 
means  to  efficiently  create  the  arbitrary  phased  sound  field  in  a  thermoacoustic  device. 
Rotating  acoustical  waves  in  cylindrical  and  ring  resonators  were  experimentally 
observed.  Single  point  excitation  was  achieved  in  both,  as  well  as  observation  of 
acoustical  shocking  in  one  of  the  resonators.  Spiral  wave  fields  were  also  explored. 
Finally,  a  thermoacoustic  cryocooler  was  built  and  tested  under  various  phases  of  the 
sound  field  impedance. 
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